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Methods of computing the mass-flow rate of nonperfect gases 
are discussed. Data for computing mass -flow rate through a 
sonic-flow nozzle are given for air, nitrogen, oxygen, normal 
and para -hydrogen, argon, helium, steam, methane, and natural 
gas. The pressure ranges to lOOxKr N/m (~ 100 atm). For 
steam, the temperature range is from 550 to 800 K. For the 
other gases, the temperature range is from 250 to 400 K. 

INTRODUCTION 

When flow meters are used for measuring the mass-flow rate of gases, 
errors may arise if the flow relations that are used in data reduction 
involve the assumption that the gas is perfect. For this report, a 
perfect gas is defined as one having an invariant specific heat and a 
compressibility factor of unity. A perfect gas is to be distinguished 
from an ideal gas , which has a temperature-dependent specific heat and 
unity compressibility factor. A nonperfect gas is a real gas . The 


assumption that the gas is perfect is sufficiently accurate for computing 
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the flow of such gases as air and nitrogen at atmospheric pressure and room 
temperature. However, for gases at high pressure or low temperatures, 
significant errors will result if the perfect-gas flow relations are used. 

There are a number of cases where the real-gas corrections are simple 
to apply. For these” cases, the changes in pressure and temperature of the 
gas as it flows through the meter are much smaller than the respective 
absolute levels of pressure and temperature; the flow can then be considered 
incompressible. It is then only necessary to know the gas density, which can 
be determined from the pressure and temperature. In references [1, 2, and 3] 
the densities or compressibility factors of some common gases are tabulated 
as functions of pressure and temperature. For these cases, the real-gas 
correction consists of using an accurate value of density rather than a 
value that would result from assuming the gas to be perfect. The correct 
density can be calculated from an equation of state, or else obtained from 
a tabulation, such as references [1, 2, or 3]> of density or compressibility 
factor as a function of pressure and temperature. Two examples where 
incompressible flow may be assumed are: 

1 . A volumetric flowmeter such as a turbine-type meter where the 
pressure drop across the meter is much smaller than the absolute 
level of pressure. 

2 . A head-type meter such as a nozzle or orifice operating at a high 
pressure level where the pressure drop across the meter is much 
smaller than the absolute level of pressure. 

The material presented in this report applies to head-type flowmeters 
through which the flow can be considered one -dimensional and isentropic. 
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Two such flowmeters are the nozzle and the venturi. In both of these meters, 
the flow from the upstream plenum to the flowmeter throat can he considered 
one-dimensional and isentropic to a good approximation. Actual deviations 
from these conditions can he handled hy applying a multiplying factor (the 
discharge coefficient) that is almost unity, and is a function of Reynolds 
number. The flow through an orifice is neither one -dimensional nor isentropic 
to a sufficient degree to permit rigorous application of the data developed 
here. 

While the conventional isentropic flow equations apply to a perfect gas, 
a number of investigators have considered the isentropic flow of nonperfect 
gases. References [4 and 5 ] develop equations for calculating the isentropic 
flow of gases described by the Van der WaaHs equation of state. A method of 
calculating the flow of nonperfect diatomic gases using Berthelot's state 
equation is described in reference [6]. In references [ 7 , 8, and 9L the 
authors consider the flow of gases described by the Beattie-Bridgeman 
equation of state. In addition, reference [ 9 ] presents tables of functions 
that aid in the one-dimensional flow calculations of real air. References 
[10 and 11] present methods of estimating isentropic exponents to be used 
for calculating the flow of imperfect gases through sonic-flow nozzles . 

(A sonic-flow nozzle is one in which the throat velocity equals the local 
speed of sound. A sonic-flow nozzle has also been called a choked nozzle 
or a critical-flow nozzle .) Reference [12], using the tabulated data of 
reference [1], presents a graphical method for computing the isentropic 
mass flow rate of imperfect gases. In reference [ 13 ], the author reviews 
the sonic -flow meter and suggests methods for correcting for gaseous 
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imperfections. In reference [14], a computer program for calculating the 
one-dimensional flow of imperfect gases is given. The program involves the 
interpolation of sets of thermodynamic-property data that are stored in the 
computer. These same authors have published a set of tables (ref. [15] ), for 
calculating the one-dimensional flow properties of real air. The thermo- 
dynamic data for air that are involved in reference [15] are the data 
tabulated in reference [1]. Sonic-flow functions for steam are given in 
reference [16] . Reference [IT] presents a set of graphs that permit the 
mass-flow rate of air, Kg, Og, n-Hg, ^ ant ^ HgO " t: ^ rou S^ L sonic-flow nozzles 
to be calculated. Information is also given on how to make the mass-flow 
rate calculation when the nozzle is operated subsonically. In reference [2], 
the air, Kg, Og, n-Hg, and HgO data of reference [17] are presented in 
tabular form; p-Hg data are also presented. Reference [3] has tabulated data 
that permit calculating the mass-flow rate of Ng and He through sonic-flow 
nozzles. Reference [3] differs from references [2 and 17] in that the 
pressure and temperature ranges in reference [3] are greater than the ranges 
in references [2 and 17]- References [2 and 3] also contain tables of such 
thermodynamic properties as compressibility factor, specific heat, and 
speed of sound. In addition, reference [3] contains the computer programs 
used in making the calculations. In reference [18], data for calculating 
the flow of natural gas through sonic-flow nozzles are presented; the 
computer program for calculating these data is given in reference [19] • 

In this report, the sonic flow data for various gases, as presented 
in references [2, 3; 17 > and 18], are summarized. Since these references 
were published, more exact calculations for argon and methane have been made. 
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The more exact argon data replaces that in reference [17], and the more 

exact methane data replaces that in reference [18] . These data are presented 

in terms of a sonic-flow factor. The use of this factor permits the mass 

flow rate of these gases through sonic-flow nozzles to he calculated. In 

addition, the empirical method, given in reference [17], of calculating the 

mass-flow rate of these gases through subsonic nozzles, is presented. For 

all gases except steam, the calculations are for temperatures from 250 to 

400 K, and pressures to 100xl0 ? w/rn . For steam, the temperatures are from 

5 2 

500 to 800 K, and the pressures to lOOxKr W/m . S.I. i960 units are used 
throughout this report. 
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SYMBOLS 

2 

area, m 

speed of sound, m/ sec 

2 2 

specific heat at constant pressure, m /(sec K) 

2 2 

specific heat at constant volume, m /(sec K) 

2 

mass flow rate per unit area, kg/(m sec) 
enthalpy, m /sec 

integration constants for enthalpy, K 

integration constants for entropy 

mass flow rate, kg/sec 
2 

pressure, N/m 

2 2 

gas constant, m /(sec K) 
entropy, m 2 /(sec 2 K) 


T 


temperature, K 
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v velocity, m/ sec 

Z compressibility factor 

T specific -heat ratio 

JL 

4> sonic-flow factor defined implicitly by equation 12, (sec K 2 )/m 

Subscripts 

0 refers to plenum station 

1 refers to nozzle throat station 

i refers to ideal-gas condition 

P refers to perfect-gas condition 


ANALYSIS 

The conditions assumed in this analysis are as follows: The gas 

is at rest in a plenum and accelerates one-dimensionally and isentropically 
to the throat of a nozzle where its speed is sonic. The measured quantities 
are the plenum pressure and temperature. The gas is not assumed to be 
perfect, and its state equation is given by 

p = ZpRT (l) 

where Z is the compressibility factor and may be expressed as a function 
of pressure and temperature or of density and temperature. The assumption 
that the flow is one-dimensional and starts from rest is represented by 
H q = H x + i v^ (2) 


The assumption that the flow is isentropic is represented by 



1 


(3) 


and the fact that the flow is sonic is represented by 


v i = a i 


(4) 



T 


In order to solve equations (l) through (4), it is necessary to express 
enthalpy, entropy, and the speed of sound in terms of either pressure and 
temperature, or density and temperature, depending on the form of the 
state equation. 

Case I. Z = Z(p,T) 

The expressions for enthalpy and entropy are integrated forms of Eqs. 
(4) and ( 5 ) in reference [IT]. 



The temperature integrals in Eqs. (5) and (6) are indefinite integrals whose 

and Kg. The values of 
and Kg depend on the choice of the gas reference state. Equations ( 5 ) 
and (6) also involve the ideal-gas specific heat. The term ideal-gas refers 
to a gas whose compressibility factor is invariant, with a value of unity; 
however, unlike a perfect gas, the specific heat is a function of temperature. 
This condition is approached as the pressure of the gas approaches zero, 
providing dissociation does not occur. The expression for the speed of sound 
is found in reference [IT]. The value of a is obtained from 

Z 2 RT/a 2 = Z - p (|) t - [z + T (§§)J/C y R ) (T) 


constants of integration are included in K^ 



where 
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If these expressions for enthalpy, entropy, and speed of sound are substituted 
in Eqs. ( 2 ), (3)> and (4), and the iteration procedures given in reference 
[IT] are then applied, solutions can "be obtained for the velocity, pressure, 
and temperature at the nozzle throat. Knowing the pressure and temperature 
at the nozzle throat, the corresponding density is determined through Eq.. (l). 
The mass flow rate per unit area through the sonic-flow nozzle then becomes 

G 1 = Pl v l (9) 

Case II. Z = Z(p,T) 

The expressions for enthalpy and entropy are given in reference [3] and 

are: 




where 
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These expressions for enthalpy, entropy, and speed of sound are substituted 

in Eqs. (2), (3), and (4). Equations (l) through (4) are then solved for 

</ 

plenum density, and for density, temperature, and velocity at the nozzle 
throat. The iteration procedures involved in this solution are found in 
reference [3]. The expression for the mass-flow rate per unit area through 
the sonic-flow nozzle is again given by equation (9). 


RESULTS AND DISCUSSION 
The Sonic-Flow Factor 

The mass flow rate of gases through sonic flow nozzles can be expressed 
in terms of a sonic-flow factor, <t>, as follows: 

m = • <t> • ( Pq/Tq ) (l^) 


For a perfect gas, the value <t>p of this sonic -flow factor is independent 
of pressure and temperature and is given by 




(r p +i)/(r p -i) 


(15) 


where T p is chosen to be 7/5 for air, nitrogen, oxygen, normal hydrogen, 
and para-hydrogen; 5/3 for argon and helium; and 4/3 for steam, methane, 


and natural gas. 


For an ideal gas, the sonic-flow factor has a value <tu that is 


a function of temperature and is given by 
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<t> 1 = lira * (16) 

P“^0 

Bae results of the real-gas calculations are presented graphically in 
Figs, (l) to (3). Instead of plotting the sonic-flow factor as a function 
of pressure and temperature the ratio of the real-gas sonic-flow factor to 
the ideal-gas sonic-flow factor is plotted. The reason for this is that, 
if the sonic-flow factor itself were plotted, the families of curves would 
cross each other for some of the gases, making the graphs difficult to 
read. In terms of the ordinates in Figs, (l) to (3), the expression for 
the sonic mass-flow rate is 



where the values of are given on the figures to which they apply. 

In the event that it is desired to use nozzles or venturies subsonically, 
but the velocities are such that the flow has to be considered compressible, 
the following equation, derived from Eq. (27) in reference [ 1 J ], applies 



Equation (18) is not based on theory, but is an approximation to actual 
subsonic calculations. For the pressure and temperature ranges involved, 
Eq. (l8) reproduces the subsonic calculations to within ^ percent for the 
gases considered in this report. The values (<K/<i>p) are given on the 
figures to which they apply. The values of the compressibility factor 
are given graphically in Figs. ( 4 ) to (6) for the gases considered in 
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this report. The perfect-gas mass-flow rate lip is given by 


( Ap/ ) 2 


2TV 


V 1 


RT„ 



1 


p (r p -D/r p ' 

V 


(19) 


L -> 

In Figs, (l) tq, ( 3 ), ( ) is plotted as a function of P and T 
for air, N p , 0^,, n-H^, p-H^, ^ ® e > HgO, ant ^ na ^ ura ^- § as> ^ ne "t* 16 

interesting results in Fig. (2a) is that even though the sonic-flow factors 
for n-Hg and p-H^ are different, the ratio (<t>/<k) is the same; further, 
over the range of pressures and temperatures considered, this ratio is 
independent of temperature. 

The natural gas data in Fig. ( 3 b) are for a particular composition. 
Therefore, this data would not apply strictly to a natural gas of a different 
composition. However, since methane is usually the principal component of 
natural gas, the sonic-flow factors of natural gases approximate those of 
methane . 

Table I lists the sources of the data presented in Figs, (l) to ( 3 ), 
as well as the references for the compressibility-factor and ideal-gas 
specific-heat data that were used in the calculations. 

The pressure and temperature ranges covered by some of the references 
exceed the ranges covered in this report. Table II lists the actual ranges 
covered in references [2, 3; 17 > and 18 ] . 


Compressibility-Factor Data 

In order to use Eq. (l8), it is necessary to have compressibility-factor 
data. To this end, Figs. (4) to (6) are presented. Pressure and temperature 
ranges are the same as in the sonic-flow ratio figures (Figs, (l) to ( 3 ))* 
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Table I. Data References 


Gas 


Sonic -Flow 
Factor 


Compressibility- 

Factor 


Ideal-Gas 
Specific Heat 


Air 

2,17 

1 

1 


2,17 

1 

1 

°2 

2,17 

1 

1 

n-H 2 

2,17 

20 

20 

p-h 2 

2 

20 

20 

A 

(a) 

21 

21 

He 

3 

22 

22 

H 2° 

2,17 

23- 

1 

ch 4 

(a) 

24 

25 

Nat. Gas 

18,19 

18 , 19,26 

1,25,27 


(a) These data are presented in this report for the first time. 
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Table II. Pressure and Temperature Ranges Covered in References. 
Pressures in 1C>5 N/m 2 (~ atm); Temperatures in K. 


Gas 

Ref. 2 
( Tables ) 

Ref. 3 
(Tables) 

Ref. IT 
( Graphs ) 

Ref. 18 
(Tables) 

Air 

Pq<100 

22°<T^390 - 

— 

Pq<100 

22°<T 0 <390 

— 


p 0 <l° 0 

220 <Tq <390 

Pq<300 

100 <T 0 <4 00 

p 0 <ioo 

22°<t 0 <390 

— - - 

°2 

Pq< 100 

22°<T C )<39 0 

— 

p^ioo 

22°<T^<390 

— 

n-H 2 

p^ioo 

100 < T o<390 

— 

Po< 100 

220<T 0 <390 

— 

p-h 2 

Pq<100 

100 < T o<39P 



— 

A 

— 

- - - 

Pq<100 

22°<To<39 0 

— 

He 

— 

Pq<300 

!5< T o<^ 00 

P 0 <1°° 

22°<T 0 <390 

— 

H 2° 

Pq<300 

4°°<t 0 <83 0 

— 

Po<3° 0 

440<To< 830 

— 

ch 4 

— 

— 

— 

P 0 <69 

2 5°< T o<390 

Nat. 

gas 

1 

— 

— 

— 

P^ 6 9 

250<T 0 <390 



FIGURE LEGENDS 


(a) Air. 

(b) Nitrogen. 

(c) Oxygen. 

Figure 1. - Ratio of the real-gas sonic-flow factor to the ideal-gas sonic- 
flow factor for various gases. 

(a) Hydrogen. 

(b) Argon. 

(c) Helium. 

Figure 2. - Ratio of the real -gas sonic-flow factor to the ideal -gas sonic-flow 
factor for various gases. 

(a) Methane. 

(b) Natural gas (fractional composition by volume, CH^-0. 9272, 02^-0.0361, 
C 3 H 8 -0.0055, iC 4 H 10 -0.0007, nC^-Cl 0010, N 2 -0. 02 18, CO 2 -O.0O77). 

Figure 3. - Ratio of the real -gas sonic-flow factor to the ideal -gas sonic - 
flow factor for various gases. 

(a) Air. 

(b) Nitrogen. 

(c) Oxygen. 

Figure 4. - Compressibility factor for various gases. 


(a) Normal and para hydrogen. 

(b) Argon. 

(c) Helium. 

(d) Steam. 

Figure 5. - Compressibility factor for various gases. 

(a) Methane. 

(b) Natural gas- (composition is the same as in fig. 3(b)), 
Figure 6, - Compressibility factor for various gases. 







SONIC -FLOW FACTOR RATIO, 9/9 j 


PLENUM 9 . 9 

j/^p 

| 

TEMPER- 1/9 

ATURE, SECK i/2 /m 

V 


PLENUM 

TEMPER- 


PLENUM 

NORMAL HYDROGEN 

PARA -HYDROGEN 


TEMPER- 

*1 

<Pj IVp 




ATURE, 

K 

SEC K 1/2 /m 


SEC K 1/2 /m 



250 

0.01073 

1.006 

0.01055 

0.989 


275 

.01071 

1.004 

.01058 

.992 


300 

.01069 

1.003 

.01060 

.994 

PLENUM 





.997 

.998 

TEMPER 

ATURE, 



(A) HYDROGEN. 







9j = 0.01593 

9 j / 9 p = 1. 000 


(B) ARGON. 


(C) HELIUM. 


PLENUM 

<P\ 


-TEMPER- 

SEC K 1/2 /m 

ATURE, 
— K 


600 _ 

/ 

550 

0.03113 

0.993 / 

— 600 

.03108 

.992 / 

650 

.03103 

.990 / 

















SONIC FLOW FACTOR RATIO, 9l<P\ 


1. 10 


PLENUM 


<p\i<pp 


TEMPER- 

ATURE, 

SEC K 1,2 /m 



K 




250 

0.02951 

0.998 


275 

.02946 

.996 

PLENUM 

300 

.02940 

.994 

TEMPER- 

350 

.02925 

.989 

ATURE, 

400 

.02907 

.983 

K 



250 275 

300 




(A) METHANE. 

NJ 

250 

.03054 

.996 

a£ 

o 

275 

.03049 

.994 

s 

300 

.03042 

.992 

£ 

350 

.03025 

.986 

£ 

400 

.03006 

.980 






o 

o 


(B) NATURAL GAS (FRACTIONAL COMPOSITION BY VOLUME, 
CH 4 -0. 9272, C 2 H 6 -0.0361, C 3 H g -0.0055, iC 4 H 10 -0.0007, 
nC^io-O-OOlO, N 2 -0.0218, CO 2 -0.0077). 


TEMPER 

ATURE, 




PRESSURE, N/m 2 


(C) OXYGEN. 





